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In the present study, a multi-criteria sizing function (MCSF) is proposed for designing the optimum size 
and operating strategy of the prime mover of a residential micro-combined cooling heating and power 
(CCHP) system. The CCHP prepares the electrical, thermal, cooling, and domestic hot water demands of 
the same building in five different climates in Iran. The MCSF integrates fuel energy saving ratio (FESR) 
and exergy efficiency as the thermodynamical parameters, net present value, internal rate of return and 
payback period for the economical criteria, and CO 2 , CO and NO* reduction for the environmental 
evaluations. Analytical hierarchy process is used to weigh each criterion with respect to others. The 
engine proposed by MCSF results in considerable fuel saving and pollution reduction and a payback 
period of about 6 years for the 5 climates. In addition, the best strategy according to the engine size is 
determined for every climate. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

A combined cooling heating and power (CCHP) system produces 
power and uses the waste heat for cooling and heating purposes. 
With respect to the conventional systems of cooling, heating and 
power generation, the CCHP systems save energy, make economical 
revenue, and reduce environmental pollution. They also have 
higher safety in crises such as war or earthquake, and flexibility to 
be on-grid or off-grid [1—6]. 

Main components of a CCHP system comprise of the prime 
mover, the heat exchangers for the heating purpose, thermally 
activated cooling system, control system and the energy resource 
used. A heat storage system may be used occasionally for storing 
the extra heat to be used in the heating or cooling energy peak 
hours. The prime mover can be a gas turbine, reciprocating internal 
combustion engine, Stirling engine, steam turbine, micro gas or 
steam turbine and fuel cells. The thermally activated cooling system 
may include technologies such as absorption chiller, adsorption 
chillers, ejector refrigeration system and desiccant dehumidifiers. A 
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compression chiller may be used in addition to the thermally 
activated systems [1,5,7—11]. 

The prime mover size of the CCHP cycle tells us what strategy is used 
to provide the electricity or heat demand of the building. In the litera¬ 
ture, three strategies are usually introduced. They are as follow [ 12—18] : 

1 Following electrical load (FEL): In this strategy, the prime 
mover nominal size (£ nom ) equals the electrical peak demand 
(Edem) and no electricity is bought from the grid. 

2 Following thermal load (FTL): In this case, the recovered heat 
(Qrec) from the CCHP cycle is the same as the peak heat demand 
(Qdem) and no auxiliary boiler is used. 

3 Following seasonal strategy (FSS) or hybrid electric thermal sys¬ 
tem (HETS): in this strategy a load ratio (LR) is calculated monthly 
and if LR > 1 the FEL is used otherwise FTL should be used. 

LR = Monthly E dem 
Monthly Qdem 

Although the numerical calculations show that, the third strat¬ 
egy is more efficient [12] but in practice, the third strategy is 
meaningless for residential buildings because it is impossible to use 
a single prime mover for both of the FEL and FTL and switch be¬ 
tween them monthly. The reason is that the size of the prime mover 
for the FEL and FTL strategies is very different for the residential 
building applications. In addition, if we try to use the third strategy, 
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Nomenclature 

1C 

Internal combustion engine 



ii 

Initial investment of components (Rial) 

AHP 

Analytic hierarchy process 

foM 

Operation and maintenance cost ($/l<W) 

amb 

Ambient 

1RR 

Internal rate of return (%) 

ATD 

Aggregated thermal demand 

i = o, i 

CSP and CCHP as superscript 

ATD 

The aggregated cooling and heating produced by 

max 

maximum 


sources other than engine recovered heat (kW) 

N 

Normalized 

b 

Boiler 

nom 

Nominal 

C 

Cooling load (kW) 

NPV 

Net present value (Rial) 

C' 

The cooling produced due to the heat of auxiliary 

P 

Power or electricity (kW) 


boiler (kW) 

PB 

Payback period (year) 

Ci, C 2 and C 3 Thermodynamical, economical, and environmental 

PHR 

Power to heat ratio 


criteria 

PVP 

Present value percentage 

CCHP 

Combined cooling heating and power 

PP 

Power plant as sub. 

cf 

Cash flow (Rial) 

Q 

Heat (kW) 

COP abc 

Coefficient of performance of the absorption chiller 

r 

Interest rate (%), and recovered as sub 

CSP 

Conventional separate production 

rec 

Recoverable 

D 

Hot water demand load (kW) 

s 

Supplied 

DHW 

Domestic hot water 

s 

Sub-criteria 

db 

Dry bulb 

sum, win Summer and winter 

dem 

Demand of building 

T 

Temperature (°C) 

E 

Electricity (kW) 

fgas 

Buying tariff of natural gas (Rial/m 3 ) 

e' 

Electricity bought from the grid (kW) 

£bp» fsp 

Buying and selling tariff of power (Rial/kW h) 

Em 

Emission 

TDC 

Tropical and dry in summer, cold in winter (Climate) 

eng 

Engine 

TDEC 

Temperate and dry in summer, extremely cold in 

er 

Earned money (Rial) 


winter (Climate) 

ex 

Expanse (Rial) 

THC 

Temperate and humid in summer, cold in winter 

EXIR 

Exergy increase ratio (%) 


(Climate) 

F 

Fuel (kW) 

THT 

Tropical and humid in summer, temperate in winter 

FA 

Fuzzy algorithm 


(Climate) 

FCU 

Fan coil unit 

TSHT 

Tropical and semi-humid in summer, temperate in 

FEL, FTL 

Following electrical and thermal load 


winter (Climate) 

FESR 

Fuel energy saving ratio (%) 

W 

Work rate (kW) 

FLO, PLO 

Full and partial load operation 

wb 

Wet bulb 

FSS 

Following seasonal load 

wh 

Water heater 

g 

grid 

wp, sp, rp Water, solution and refrigerant pumps 

GIA 

Grey incidence approach 

y 

Yearly 

H 

Heating load (kW) 

VCHP 

Combine heat and power efficiency 

H' 

The heating produced by the auxiliary boiler (kW) 

Coverall 

Overall efficiency (%) 

hg 

Hot gas 

7r 

Exergy efficiency (%) 

i 

Investment or emission index (Rial/kW h), (Rial/m 2 ), 

0 

Exergy rate (kW) 


(kg/MW h) 

A, V 

AND, OR logical operators 

i 

Total investment (Rial) 

0 

Thermodynamical dead state as sub 


the prime mover must operate in partial load some months, which 
results in efficiency reduction. 

After selection the type of prime mover, sizing of prime mover is 
the most important parameter. Sizing changes the amount of fuel 
saving, overall efficiency, environmental pollution, economical 
parameters and the optimum strategy used to prepare the electrical 
or thermal loads of building. 

Ref. [19] uses two multi-criteria decision-making methods to 
select the best type of prime mover for five climates in Iran. The 
authors use 16th sub-criteria from the technological, economical, 
environmental and social criteria for the comprehensive decision¬ 
making. The effect of climate on sizing and designing of CCHP is 
studied by Ref. [6], The authors used maximum rectangle method 
(MRM) to determine the prime mover size of the CCHP system for a 
similar residential building in five different climates of Iran. 

Ref. [20] determines the size and number of micro-gas turbines 
of a CCHP system for a residential building in Tehran, Ahwaz and 
Hamedan when using the FTL strategy. This research uses a simple 
and linear balance equation between the peak loads of the building. 


A thermo-economic analysis is presented by Ref. [14] to specify the 
type and number of required micro-gas turbines to provide the 
electrical and thermal needs of a commercial building in Tehran 
urban area during a year. This research maximizes an objective 
function called annual profit. Refs. [18,21] have used the MRM for 
sizing the components of the CHP or CCHP systems. The principals 
of MRM are presented in Refs. [14,18,21]. Ref. [21] considers a 
constraint on the primary energy saving as the main criteria and 
using FTL strategy. They used the hourly cooling load for sizing the 
absorption chiller and the aggregated thermal demand to size the 
prime mover of a CCHP plant. In 2010, Ref. [18] developed the MRM 
presented in Ref. [21] to three methods of vertical design, hori¬ 
zontal design and high-level analysis for sizing the prime mover. 
The only constraint in the vertical and horizontal design is a min¬ 
imum pre-assumed PES, but in the high-level analysis an 
economical criteria such as net present value (NPV) is also 
considered. 

In the present paper, a novel method that considers energy, 
exergy, economical and environmental aspects of the project is 
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presented for sizing and strategy selection of the prime mover of 
the CCHP system. By linear combination of the criteria and calcu¬ 
lating weight for each criterion a function called multi-criteria 
sizing function (MCSF) is achieved. This function is used for opti¬ 
mum sizing of the prime mover from the single or multi criteria 
points of views. The results of sizing are presented for five different 
climates by considering a similar and hypothetical but very close to 
real residential buildings in Iran. 

2. Formulation 

2.3. Calculate the criteria and sub-criteria 

In this research, thermodynamical (Ci), economical (C 2 ) and 
environmental (C 3 ) criteria are considered for finding the best prime 
mover size and strategy for the CCHP system in different climates. 
Exergy and energy analyses are done for the thermodynamical an¬ 
alyses. The economical evaluation is performed by using NPV, in¬ 
ternal rate of return (1RR) and payback period (PB) to find out about 
the profitability, risk and payback period of investment in CCHP. For 
the environmental assessment of the CCHP cycle, productions of the 
NO*, CO 2 and CO are considered. In the following, the equations used 
for different evaluations are derived systematically. The magnitudes 
of the criteria are calculated for both the CCHP cycle and conven¬ 
tional separate production of heating, cooling and power (CSP). The 
cooling system in both CCHP and CSP is an absorption chiller. An 
auxiliary boiler and recovered heat from the CCHP cycle are used for 
providing the heating load of the building and domestic hot water. 
Fan coil units (FCU) are used to transfer heat and cool to the final 
users. The schematic of the CCHP and CSP are shown in Fig. 1. 


The components of fuel vector for the CSP and CCHP are pre¬ 
sented below: 


pO _ Qlem pi _ E + Qri 


VppVg 

pO _ Qlem 
C tibCOPabcfipcu 


4CHP 


U 1 _ 

— 




ec | (^dem) 


VppVg 


"dem 


^bCOPabcepcu 


pO _ ^dem pi _ 1 Mem 
H Vb £ ¥ai H ^b^FCU 

cO _ ^dem cl _ Adem 
D e wh ’ ° ^wh 

In which: 


(4) 


E = 

Qrec — 
Qrec = 

F' = 

L dem 

C' = 

dem 

n dem “ 


£nom if(£nom < ^dem) V FLO V [(£ n om < ^demJ^PLO] 
^dem if[(E nom > £dem)^PLO] 

1.368E + 14.57,30 < E(kW) < 400 
1.854E, 0 < F(kW) < 30 
J Qlem — Pnom if Pdem > Pnom 
( 0 if P,] ern < Pnom 

Qlem — COPabcQiec if Qlem > COPabcQrec 
0 if Qlem - COP d |, c Qjec 

Hdem ~ Qrec if FQ em > Qrec 
0 if Pf,] em < Qrec 


(5) 


2.3.3. Thermodynamical analyses 

Thermodynamical analyses include the energy and exergy an¬ 
alyses. In the energy analysis, the fuel consumption of the CCHP is 
compared with the CSP. However, in the exergy analysis the overall 
exergy efficiency of the CCHP is compared with the overall exergy 
efficiency of the CSP. 

In the energy analysis, a fuel vector (F) is defined as below: 

F> = (/"p- P( • I'h■ P|>) (2) 

in which the subscripts of P, C, H and D stands for the power, 
cooling, heating and domestic hot water and superscript) = 0 and 1 
are for the CSP and CCHP respectively. To calculate the fuel con¬ 
sumption during a year we have 8760 h in one year, hence: 


Attention must be paid to the water heating system for domestic 
hot water (DHW) demand. If it is supposed to use the recovered 
heat for preparing the DHW, therefore Ddem = 0 and 
Hdem = (heating load + DHW load). However, if a separate water 
heater is installed Ddem = DHW load and Hd em = heating load. In 
addition, the aggregated thermal demand (ATD) of the building is 
calculated by ATD = Hdem + Qem/COP a b C , and correspondingly 
ATD' = H' dem + Cd em /COP abc shows the aggregated thermal lack of 
the CCHP system that must be compensated by an auxiliary boiler. 

In the exergy analysis, the recovered and supplied exergy vec¬ 
tors for the CSP and CCHP are calculated, and used to compute the 
exergy efficiency during a year as below: 

= ^Vs,P) $s,O 0s,Ht 0s, dJ (6) 


8760 

Fuel J = Z (4 + 4 + ^h+^d), 



(7) 



Fuel+Oxy gen 


Conventional 
Power plant 


Return water 

Stack | 

E a „. 1 r 


dem 
^ ^ Lost 


Boiler 
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Fuel+Oxygen 

Suck | Return water Supplied Chilled water 


■ Supplied I 


Absorption chiller 


Supplied Hot water 
Fuel+Oxygen 

W nilmt 


Ireul - 


Fig. 1. The CCHP system (left) and CSP system (right). 
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nj = E /£ (V*) (8) 

t=i V /t t= i V /t 

In which fe is the number of time steps in one year used for the 
load calculation, In addition, the components of the exergy vectors 
are calculated according to the following equations: 

4>s,c = Wsp + vv rp + W wp + F° - F°(l - Vb ) (l - 

4-vv„ P+fS ( 

(9A) 


j amb 


R db ,sum if(H d 

em — OACdem ^ 0) 
Tdb,win if(Qiem = 

Troom if f^dem Qjem = ^ 


(1 IB) 


Troom.sum 

Troom,win 

(7room,sum + Troom,win)/2 


if(Qlem — 

if(ffdem = O^Qlem 9 ^) 
if H dem = C dem = 0 

(11C) 


The input data for the thermodynamical analyses are given as 
follow; r/pp — 0.3, ijg — 0.85, ??b — 0.8, epcu - 0.85, eyvh — 0.85, 
^/CHP = 0.8, COPabc = 0.7 VV S p + l^rp = 0.007C n om(kW), W wa ter pump 
= 0.005max(C dem , H dem )(kW), T 0 = 15°C, T hg = 1000 oo C, 
^oil — 00 °C, Tjacketing — 100°°C, T ex haust — 540°°C, T r0 om,sum — 25°°C, 
Troom.win — 22°°C, T d b,sum, Tehran — 38°°C, T d b,win, Tehran — — 2.9°°C, 
Tdb.sum, Ahwaz = 46.74°°C, T d b,win, Ahwaz = 8.73°°C, T d b,sum, Bandar Anzali 
— 29.76°°C, T d b,win, Bandar Anzali = 4.61 °°C, T d b,sum, Chabahar = 33.55°°C, 
Tdb.win, Chabahar — 16.87°°C, Tdb.sum, Kmayaran — 36.82°°C, Td^win, 
Kamyaran — —4.88°°C. 


rr,P 


= E, 


dem 


rr,C 

,0 


= c 


dem 


1 


r.H = W dem (l- T ° ) 
1 room / 


2.1.2. Economical analyses 

In the economical analyses, the NPV, IRR and PB of the CCHP 
system are calculated. The initial investment is a key parameter in 
the computation of the NPV and IRR. The following equation is 
(9B) recommended for calculation of the initial investment cost (I) for 
the CSP and CCHP systems: 



P = iP p + i*c + iib + “d ( 12 ) 

In which the ii is the initial investment cost for each part of the 
CSP or CCHP systems. 

iip = Subscription fee to grid, iip = E nom -ip (13) 


4>s,c Wsp + lV rp + VV w p + Fq ^1 — — R c(l — Pb) ^1 ~ j ~) 

K„ = H'«P+Fi(l-^)-F|5(l-1b)(l- i 5 ; ) 

(10A) 


ii? — ^nom'^C’^C — f-nom ’Iq 
Qiom = max(C dem ) 

ilH — TTnom'^H’^H = (Tfnom — Qrec) ‘ I'h 
TTnom max(TT de m? f'dem/^'f^f > abc) 

iip = Dnom*^D)iiD = E) nom ij^ 

Dnom — max(D d em) 


(14) 

(15) 

(16) 


0r p - E + E dern + 00^ 

\ iroom/ 


r 'Vl I Ale m 


1 - 


To 


K, D = D d em( 1 - t T ° ) 
*room/ 


(10B) 


In which, the recovered heat sources include the exhaust gases, 
water-jacketing and oil cooling system, it is assumed that they 
approximately account for 50%, 35% and 15% of the total recovered 
heat from the engine respectively. In addition, T am b is defined for 
each city according to the weather data [19,24]. 

Qrec = Qoil ‘ Qjacketing ' Qexhaust 

^ = Q °" ( 1 “ 3 + QjaCketinS “ Acting) ( 11A ) 

T Qexhaust (1 — r 

V ' exhaust 


In which i is the index price (Rials/kW) for each part of system. 
To calculate the annual cash flow (cf) for the CSP and CCHP systems, 
the annual expanse vector (exv) is subtracted from the annual 
earning vector (erv) as below: 

eiV = (Wp,er J c ,er J H ,er J D ) (17) 

exV = (exj,, ex J c , ex J H , ex J D J (18) 


8760 

cf J y = ^2 (^erV - exV j 


n 


cf!, 


NPV J = V- y -^ - P 

y A1 (1 + ry 


n c fj 

y' y 
(i + IRR ) y 


= p 


(19) 


( 20 ) 

( 21 ) 
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In which, n — 30 is the lifetime of the project. Also, the erv 7 and 
exv 7 are calculated as below: 


er p — 0 ? — E dem ' ^bp 

er° = 0, exg = F° ■ t gas 
er H = 0,ex° = Fg-tgas 

er D = 0ex S = f D f gas 


( 22 ) 


er 

er 

er 

er 


i _ 

p — 


i _ 
c 


i _ 
H 


0 _ 
D — 


Fdem ' fbp + (F F dem )' Fp ■ 6xJ 

Qr ec . f PX 1 _ pl.t- 
—— l gas ■ eA c — r c l gas 


wx.t PX 1 _ F 1 

—— tgas,ex H — r H L gas 

'lb 

0 ex f u = F’-tgas 


Fp •tgas 


(23) 


The third criterion that is usually used in the economical ana¬ 
lyses is the payback period of the project. This criterion calculates 
the years that it takes to recover the original initial investment cost, 
and is calculated as below: 


PB' = P/cPy 


f PB' if cf y > 0 

\ not applied if cfy < 0 


(24) 


The input data for the economical analyses are given as follow; 
t bp = 485 Rials/kW h, t sp = 392 Rials/kW h, 

tgas—(1 - offccHp) x 700 Rials/m 3 , offccHP = 20%, iiJS = 
3.45 x 10 6 Rials, in which offccHp is the government discount for 
selling natural gas to the CCHP users. In addition, the equations for 
calculating the components’ price are presented in Table 4 [25]. 


2.1.3. Environmental analyses 

In the environmental analyses, an environmental vector is 
defined for the CO 2 , CO and NO* production in the CSP and CCHP 
systems as follow: 


2.2. Normalization of the criteria 


In order to normalize the sub-criteria of S, if it is a cost or “the 
lower the better” like the fuel consumption, then the reduction 
ratio of S can be calculated as below: 

= lOO^-S^/s 0 (28) 

However, if S is a benefit or “the higher the better”, then the 
increase ratio of S can be calculated as below: 


= lOO^S 1 -S 0 )^ 1 (29) 

Eqs. (28) and (29) are used to normalize the technological and 
environmental sub-criteria. In order to normalize the economical 
sub-criteria of NPV, IRR and PB the following techniques are used: 

The NPV of separate production is always negative, because this 
system never has positive cash flow for the owner; also only the 
CCHP projects with positive NPV is economically acceptable. Hence, 
the NPV increase ratio is calculated as below: 


NPV" = 100- NPV 


NPV 1 - NPV 0 


f NPV" if NPV 1 > 0 
1 0 if NPV 1 < 0 


(30) 


For an economical project, the IRR shows the risk of investment, 
and must be greater than the interest rate (r). Due to the negative 
NPV 0 , the IRR 0 is always smaller than the r, for this reason the IRR 
increase ratio is calculated as below: 


IRR" 


100 


IRR 1 - r 
IRR 1 


f IRR" if IRR 1 > r 
1 0 if IRR 1 < r 


(31) 


The third economical criterion, PB, cannot be defined for the 
separate production system because there is no positive cash flow 
for the separate production. Therefore, it is normalized as below to 
be used in the MCSF: 


PR" = 122 = J PB" if PB 1 > 1 
PB 1 \ 1 if PB 1 < 1 


(32) 


x j = (x|,x'..x; l .x',j 


(25) 2.3. Conditions and constraints to find £j{ om 


In which X can be CO 2 , CO or NO*. In addition, the components of 
the environmental vector are as below: 


Xp — F dem ■ 'x.P’Xp — Bi xp + Extern’ 

'x,c 


X° = C, 


dem''x,C’X(; — C dem ! '’ 


‘X,P 


xg = H dem iUXu = H[ 


X° = D, 


dem'x.D' d 


,xA = 


;i 

'X,H 
B'dem' 'x,D 


(26) 


dem 


Therefore, the environmental emission production for each 
pollutant (Em) during a year for CSP and CCHP systems can be 
calculated according to the following equation: 


In order to generate the MCSF, we need to calculate the engine 
size (Fnom) recommended by each criterion. For this purpose, the 
following procedure is used to calculate corresponding E nom for 
each criterion. For the thermodynamical and environmental sub¬ 
criteria: 

Fnom = |F n om|S"iS max, i = 1,3, £mi n < E n om < FmaxJ - 

(33A) 

In addition, for the economical sub-criteria, the engine size is 
chosen for the conditions in which the project becomes profitable, 
safe, and the original investment recovers less than 6 years. Hence: 


8760 

Em/ X = E (X 1 , + + Xg + X J D ) t (27) 

The input data for the environmental analyses are given as follow 
[26,27]; i° op = 4, i° oc = 0.1274, i° 0H = 0.1274, i° 0D = 0.1274, 
'co 2 ,P = 700, i° 02C = 182.04/COP abc , i° 02-H = 182.04, f COiD = 

182.04, i° a P = 300, i° a c = 0.1532, i° a H = 0.1532, i° a D = 

0.1532,i'A op = 0.8,i'A oc = 0.1274,i'A OH = 0.1274,iJ. 0D = 0.1274, 
'co 2 ,p = ! co 2 ,c = 182.04/COP abc , ! co 2 ,h = 182.04, ! co 2 ,d = 

182.04, iJjo x ,p = 0-7, ’JjOj.c = 0.1532, iJjo x ,H = 0.1532, iJio^d = 

0.1532, the units of the above parameters is kg/MW h. 


Fnom — I Fnom 15" >0,1 — 2 ,J — l,2,F m j n < F n om < F m ax | 

(33B) 

Fnom = | F„om |S" > g , 1 = 2 ,j = 3,F m j n < F n om < Fmaxj' 

(33C) 

If a criterion (S" n ) continuously increases, as the £ n0 m increases, 
then the upper bound, £ ma x = max{£[( om , iJ=m,j=tn}. It means 
that £ max will be limited to the biggest size recommended by 
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other criteria. This constraint helps to avoid over sizing of 
engine. In addition, the lower bound is limited by the 
economical criteria to make the project economical. Then, £ min < 
minjEnom- ! = 3 ,j = 1,2,3}. 


2.4. Calculate the normalized weight matrix 


The degree of importance of each criterion (or sub-criterion) 
with respect to the other criteria (or sub-criteria) is considered by 
calculating the weight of each criterion (or sub-criterion). For this 
purpose, the pairwise comparison matrix (PC) between the n 
criteria is made as below: 


PC = [ pc d„xn (34) 

In which PCy represents the importance of ith criterion with 
respect to jth criterion. In the fuzzy algorithm, a linguistic pairwise 
comparison presented by Ref. [22] is used to provide the PC matrix 
(Table 1). However, in the Grey incidence approach (CIA) an 
equation is recommended as below: 

PCy = 1 ±0.15,5 = 0,1,2,...,10 (35) 

In which the plus sign is for the time when the ith criterion is 
5 degree more important than the jth criterion, but the minus sign 
is for the time when the jth criterion is 5 degree more important 
than the ith criterion. In this study, to simplify the pairwise com¬ 
parisons a combination of the methods used in the FA and GIA is 
used. Table 1 shows the combined method for pairwise comparison 
[19], The pairwise comparisons of the main criteria and sub-criteria 
are presented in Table 2. In this table, , S^ 2 , S 21 , S 22 , S 23 , S 31 , 
S^ 2 and S^ 3 are symbols for fuel energy saving ratio (FESR), exergy 
increase ratio (EXIR), NPV increase ratio, IRR increase ratio, PB in¬ 
verse, NO* reduction ratio, CO 2 reduction ratio, and CO reduction 
ratio respectively. The subjective weighting method of analytic 
hierarchy process (AHP) is used that is shortly presented below [4]: 


pj = E PQ/ 

i = 1 

Pmax = max (Pi,P2> ■■■;Pn) 

Pmin = rn' n (Pl 1 P2i •••! Pn) 

@ = Pmax / Pmin 

f 1 + (Pj — Pk) / Pmim 
^ik — r (3 + (Pi — Pfc)/Pmin) > 

lo 

Cik = log b ik , i, k = 1,2,..., n 
1 n 

dik = ~ ^ ~2( C H ~ c kl ) 


(36) 

Pi > Pk,0* 1 

Pj<Pk,0* 1 (37) 

e = 1 

(38) 

(39) 


Table 1 

Combine pairwise comparison of fuzzy and GIA methods. 


Fuzzy linguistic comparison of the ith criteria 
with respect to the jth criteria 

GIA 

PC 

PC(R) 

Just equal (JE) 

1 

1 

Equally important (EI) 

1.2 

0.8 

Weakly more important (WMI) 

1.4 

0.6 

Strongly more important (SMI) 

1.6 

0.4 

Very Strongly more important (VSMI) 

1.8 

0.2 

Absolutely more important (AMI) 

1.9 

0.1 


Table 2 

Pairwise comparison for the main criteria and sub-criteria. 



c, 

C 2 

C 3 



Sfi 

S" 

C, 

JE 

EI 

SMI(R) 


Sn 

JE 

EI(R) 

C 2 

E1(R) 

JE 

SMI(R) 


SN 

EI 

JE 

Cj 

SMI 

SMI 

JE 






s« 

s« 

c N 

J 23 



cN 

3 32 

c N 

3 33 

S 21 

JE 

SMI(R) 

VSMI(R) 

S« 

JE 

EI 

EI 

cN 

^22 

SMI 

JE 

SMI(R) 

cN 

d 32 

EI(R) 

JE 

EI 

c N 
d 23 

VSMI 

SMI 

JE 

c N 

J 33 

EI(R) 

EI(R) 

JE 


b\ k = 10 d * 


(40) 


Therefore, the eigenvector of the AHP is given as below: 


Ei = 


N 


n 


II *4 


i = 1,2, ...,n 


(41) 


Finally, the normalized weight matrix would be as follow: 


w,- 


Ei 

£"=t Ei 


(42) 


2.5. Calculate the MCSF 


In the present study, a MCSF is defined to optimize the size of 
the prime mover and the strategy used in the CCHP cycle. The MCSF 
uses technological (Ci), economical (C 2 ) and environmental (C 3 ) 
criteria. Each criterion consists of some sub-criteria. The MCSF is 
stated as below: 

m m 

MCSF = ,J>,- = 1 (43) 

i = l i = 1 

In which the w; and m are the weight of each criterion and 
number of main criteria respectively, and j 8 e ( 0 , 1 } is the criterion 
coefficient and is used for the purpose of single or multi criteria 
evaluations. When a criterion is supposed to be counted in the 
evaluations (3 = 1, otherwise /? = 0. The Q is also stated based on the 
corresponding E[{ om of sub-criteria of each main criterion as below: 

q = E w r £ n°m (44) 

j=l 

In which, n is the number of sub-criteria defined for each main 
criterion and, wy is the weight of each sub-criterion. The optimum 
size of the engine is equal to the magnitude of MCSF. Hence, 

E°P“ = MCSF (45) 

2.6. Strategy classification 

As discussed, the literature introduces three strategies of FEL, 
FTL and FSS. In all of the researches concerning about the strategy, 
they assume one of the strategies from the beginning and calculate 
the corresponding engine size, or they consider three types of 
strategies and according to the energy, economical and or envi¬ 
ronmental evaluations they finally choose one of the strategies as 
the best one [12—18,20,21]. In fact, they specify the size of engine 
according to the strategy assumed and the electrical and thermal 
loads of the building. It means that only three engine sizes are 
calculated for the three strategies and the best size and strategy is 
chosen among the three sizes. However, is there any guarantee that 
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the best size is confidently one of the three sizes? According to this 
research, the answer is, no. 

The authors of the present research believe that strategy is not 
an independent variable and should not be treated like that; 
because it can lead to non-optimized solutions. Strategy is 
completely dependent on the size of the prime mover, and the 
operation type (partial or full load) of the engine. The MCSF helps 
us to size the prime mover and decide about the operation type. 
Therefore, the strategy is not limited to the three types. As a 
guideline, Table 3 presents a classification for five strategy types for 
the time when the engine is operating at full load. 

3. Case study 

The residential building and the climates, which are considered 
in this study, are the same as that of Refs. [6,19], The building is 
hypothetical, has four floors, eight units, and total living area of 
1200 m 2 for 32 inhabitants. The building is considered for five 
nominated cities of Ahwaz, Tehran, Bandar Anzali, Chabahar and 
Kamyaran located in the climates of TSHT, TDC, THC, THT and TDEC 
respectively. To calculate the electrical load of building, the power 
consumptions of the refrigerant pump and solution pump of the 
absorption chiller and the water recirculation pump that were not 
considered in the Ref. [19] are considered in this research. For this 
reason, the electrical load reported in this study is a little bit 
different with respect to that presented in Ref. [19]. 

4. Results 

4.3. Load results 

In addition, according to the findings of [6] it is considered that 
the prime mover is always operating at full load. Fig. 2 presents the 
maximum ATD, maximum F^em and power to heat ratio of the 
building (PHRdem). as it can be seen due to the small changes in the 
electrical load; the PHRdem of building has changed slightly (with 
respect to Fig. 10 of Ref. [19]). The power to heat ratio of the 
building (PHRdem) varies from 0.16 to 0.34 for the nominated cities 
from the five climates; it means that the heat demand is 2.94—6.25 
times the power demand of the building. The power to heat ratio of 
the internal combustion engine (PHReng) that is considered as the 
prime mover for the CCHP system varies from 0.5 to 0.74 [23], 
which means that the Q, ec is only 1.35—2 times the £ n om- This 
comparison tells us if we follow the electrical load to provide the 
heat demand a huge auxiliary boiler is necessary, this may 
decrease the overall efficiency and exergy efficiency, and increase 
the environmental pollution, but if we follow the thermal load an 
enormous prime mover should be used that increases the initial 
investment cost of the CCHP. According to this discussion, the 
optimum size and strategy should be chosen based on the energy, 
exergy, economical and environmental analyses as follow. 


Table 3 

Strategy classification for a prime mover at full load operation. 


Strategies 

Conditions 

Esell 

Epurchase 

Lower FEL (LFEL) 

^nom < ' 

Edem.max 

0 

yi(Edem — E nom ) 

FEL 

^nom ~ ^dem.max = 

Efel 

y~l(Enom — Edem) 

0 

Following electrical, 
partial thermal 
load (FEPTL) 

Efel < Edem,max 

< Eftl 

y^fEnom — Edem) 

0 

FTL 

Qrec = Qdem.max 
Enom = Eftl = 
/(Qrec) 

y~l(Enom — Edem) 

0 

Upper FTL(UFTL) 

Enom > EpTL 

y .(Enom — Edem) 

0 


4.2. Weighting results 

The weight of each criterion with respect to other criteria is 
calculated by using the AHP method and is reported in Table 5. 
According to this table the economical criteria receives the biggest 
weight (0.434), then thermodynamical criteria (0.366) and finally 
the environmental criteria (0.200). Among the thermodynamical 
criteria, the energy weight (0.550) is slightly greater than exergy 
(0.450). Among the economical criteria the NPV weight (0.522) is 
the biggest, because this criterion talks about the profitability of 
the project while considering the interest rate (in this study 
r = 15%). The IRR weight (0.303) receives the second priority, the 
higher the IRR (with respect to the r), the smaller the risk of in¬ 
vestment. About the environmental sub-criteria, CO weight 
(0.380) is of greatest importance, because producing this toxic gas 
in the residential applications is very dangerous. After CO, the CO 2 
weight (0.331) is placed in the second order, because of its great 
impact on the global warming. The NO* weight is 0.289 and is 
placed at the third stage. 

4.3. Single criterion sizing results 

The single criteria sizing refers to finding the E[{ om for each 
criterion according to the conditions and constraints defined in the 
previous sections. The results of single criteria sizing are presented 
in Figs. 3—10 and Table 6. Fig. 3 presents the impact of changing the 
engine size on the FESR; it shows that there is a maximum of FESR 
for each climate. It shows that FESR starts increasing rapidly by 
increasing the £ nom until it reaches a maximum point and then 
decreases slowly. The reason is that as the £ n0 m increases the Qr ec 
increases as well, this increases the FESR until the Qrec reaches the 
ATD. When the Qre C becomes greater than the ATD, the heat is 
wasted and the FESR decreases accordingly. This may happen 
during days and nights, but the yearly data are presented in the 
results. The £ nom that maximizes the yearly FESR is presented for 
each climate in Table 6. Fig. 4 shows the EXIR versus £ n om. as it can 
be seen the EXIR increases as the £ nom increases. It is clear that it is 
impossible to find an £„om to maximize the EXIR. To solve this 
problem in the MCSF, the biggest £ nom recommended for other 
criteria is chosen for the criterion. For example, according 
to Table 6, the NO* reduction ratio for all of the climates recom¬ 
mends the biggest £ nom among the other criteria; therefore, the 
Enom for the EXIR is the same as £ n0 m for the NO* reduction ratio. 
This also avoids unnecessary over sizing of engine. 

Figs. 5—7 present the normalized economical sub-criteria versus 
the £ nom - The constraints for the economical evaluations are pre¬ 
sented in Table 6 as well. Positive NPV means that the project is 
profitable, IRR shows the risk of investment, and must be greater 
than the interest rate (r) as a minimum expectation. In addition, the 
maximum PB defined for this study is assumed 6 years. These 
constraints provide the minimum economical expectation of the 
designer. It is clear that these constraints can be considered more 
conservatively. For example it can be considered that NPV > 0.53o, 
IRR > (r + 5%) and PB < 5years. Figs. 5 and 6 present the NPV n and 
IRR n . Negative NPV n and IRR N are set equal to zero and only the 
positive values are economical. These two criteria recommend 
minimum engine sizes of 63,165, 66,148 and 59 (kW) for the cities 
of Tehran, Ahwaz, Bandar Anzali, Chabahar, and Kamyaran 
respectively. These sizes guarantee the profitability and safety of 
investment in CCHP. The third economical criterion concerns about 
the period it takes to recover the initial investment (lo) which is 
different for each climate. The PB N is shown in Fig. 7 and the 100/6 
line is drawn to illustrate the limit for the sizes that result in PB of 
less than 6 years. The £ nom s with PB N of greater than 100/6 mean 
payback period less than 6 years. The £ mm s recommended by PB 
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Table 4 

Cost curves of components. 

Component Equation 

Steel boiler y = 1.524 x 10 13 * 3 212 + 1.033 x 10 5 R sq uare = 0.9951 

Gas internal combustion y = 1.1 x (1.881 x 10 21 x -1033 + 1.062 x 10 7 ) +yoM, 

engine y 0 M = 0.40458 + (324,890/x)R square = 0.9562 


Absorption chiller 


y = 12,260( —35.41og(x) + 431) 


Description 

x(kW), y(Rials/kW), x < 3486, Fire tube boiler, pressure 4 bar, 
maximum temperature of 120 °C, for every 1 bar pressure 
increase, the cost increases 5%. http://www.spac.ir. 
x (kW), y (Rials/kW), x < 500,Costs are expressed in Rials per 
kilowatt of installed capacity. Data from California energy 
commission, 2004, CERTS Program Office Lawrence Berkeley 
National Laboratory, 1USD = 12,260 Rials,the 1.1 is the 
coefficient to consider the heat exchanger costs, 
x (kW), y (Rials/kW), Data from California energy commission, 
2004, CERTS Program Office Lawrence Berkeley National 
Laboratory, 1USD = 12,260Rials 



Tehran Ahwaz Bandar Chabahar Kamyaran 
Anzali 

■ PHR ■ Max. of Aggregated thermal demand ■ Max of electrical demand(Emax) 


Fig. 2. PHR, maximum of ATD and electrical demand of building. 


are 82, 217, 86, 206 and 77 (kW) for the cities of Tehran, Ahwaz, 
Bandar Anzali, Chabahar, and Kamyaran respectively. 

Figs. 8—10 represent the pollution reduction ratio of CO 2 , CO and 
NO* respectively. The CO 2 reduction ratio recommends E n0 mS of 58, 
115, 58, 93 and 57 (kW). The behavior of pollution reduction ratio is 
the same as FESR, because the pollutant production is directly 
proportional to the fuel consumption. As it can be seen the CO and 
NO* reduction ratio changes very little after a certain size of £ 110 m. 
but however a maximum point can be found for these criteria. The 
Fnom recommended by the CO and NO* reduction ratio are shown in 
Table 6. As it can be seen, the £ n0 mS recommended by different 
criteria is different and it is very difficult to use one of them as a 
reference to find the optimum size for the engine. Therefore, a 
solution should be proposed to recommend an engine size to use all 
of criteria and their corresponding engine sizes. The MCSF is the 
proposal and the results are presented in the following. 


give different E nom s. To solve the diversity of E nom s from different 
criteria and simultaneously make use of the advantages of all 
criteria a MCSF is proposed. The MCSF is calculated according to the 
weights of criteria, weights of sub-criteria, and £ n0 mS by using Eqs. 
(43)—(45). The results of multi criteria sizing are presented in 
Table 7. This table reports very important results including, 
£™m mum , FESR, EXIR, NPV, 1RR, PB, CO reduction ratio, C0 2 reduction 
ratio, NO* reduction ratio, I 0 , PVP (present value percentage) and 
the auxiliary boiler capacity for the CCHP system in the five cli¬ 
mates. The £nSm mum for different cities of Tehran, Ahwaz, Bandar 
Anzali, Chabahar and Kamyaran are 83, 211, 88,161, and 78 (kW). 
These engine sizes give annual saving energy of 36.23,35.27, 34.69, 
40.08 and 37.38% for the 5 cities mentioned above respectively. The 
EXIR for the five climates ranges from 90% to 95%. About this 
parameter, as the engine size increases the EXIR increases as well. 
When studying Table 7, it can be seen that due to considering non¬ 
conservative economical constraints, the NPV after 30 years is not 



0 50 100 150 200 250 500 350 400 

Enom 


4.4. Multi-criteria sizing results 

As discussed in the previous section, it is very difficult to use a 
single criterion to size the prime mover, because different criteria 


Table 5 

Weights of criteria according to the AHP. 


Main criteria Sub-criteria 


Technological 

0.366 

Energy 

0.550 



Exergy 

0.450 

Economical 

0.434 

NPV 

0.522 



IRR 

0.303 



PB 

0.175 

Environmental 

0.200 

NO x 

0.289 



C0 2 

0.331 



CO 

0.380 


Tehran Ahwaz Bandar Anzali Chabahar — — Kamyaran 
Fig. 3. Impact of engine size (E„ om ) on the annual S)) for 5 climates. 



Enom(kW) 

Tehran Ahwaz Bandar Anzali Chabahar — — Kamyaran 


Fig. 4. Impact of engine size (E nom ) on the annual S^for 5 climates. 
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0 50 100 150 200 250 300 350 400 

Enom(kW) 


- =Tehran - Ahwaz - Bandar Anzali - “Chabahar — — Kamyaran 

Fig. 5. Impact of engine size (E nom ) on the annual for 5 climates. 



0 50 100 150 200 250 300 350 400 

Enom(kW) 


- “Tehran - Ahwaz - Bandar Anzali —— -Chabahar-Kamyaran 

Fig. 6. Impact of engine size (E nom ) on the annual S) 2 for 5 climates. 

considerable (from 7.12MT for Chabahar to 30.04MT for Ahwaz); 
the risk margin (IRR-r) is very narrow (from 0.44% for Chabahar to 
1.67% for Bandar Anzali) and the PB is about 6 years. The 
economical constraints depend on the designer decision, they can 
be considered more conservative as discussed previously, and as a 
result, the economical criteria will improve accordingly. The envi¬ 
ronmental analyses show very considerable reduction in the pol¬ 
lutants. Green house gas, CO 2 reduces annually from 43.51% for 
Ahwaz to 46.98 for Chabahar. CO and NO* reduction of about 80 and 
90% are achieved for all of the climates. 

The initial investment, I 0 of the CCHP for the 5 climates is very 
different. It varies from 137.27MT for Kamyaran to 346.97MT for 
Ahwaz. By presenting an evaluation between NPV with the corre¬ 
sponding 1 0 , the present value percentage (PVP) is calculated. Ac¬ 
cording to the PVP that is presented in Table 7, using the CCHP in 



Enom(kW) 

Tehran Ahwaz Bandar Anzali Chabahar — — Kamyaran 


Fig. 8. Impact of engine size (£ n0 m) on the annual .S')) for 5 climates. 



0 100 200 300 400 


Enom(kW) 

- Tehran - Ahwaz - Bandar Anzali — “Chabahar ’Kamyaran 

Fig. 9. Impact of engine size (£ nom ) on the annual S) 2 for 5 climates. 

Bandar Anzali is the most economical, but in Chabahar it is the least 
economical. The CCHP in Tehran and Kamyaran is also more 
economical with respect to Ahwaz. 

The auxiliary boiler capacity for Ahwaz is very big; this is due to 
very high cooling load of Ahwaz. This big size of boiler avoids over 
sizing the engine. Figs. 11 and 12 present the effect of £ nom on the 
auxiliary boiler capacity and I 0 for the 5 climates. Fig. 11 also shows 
the engine size that result in not using an auxiliary boiler in the 
system. 

As an overall discussion on the results of this section, the initial 
investment cost for two cities of Chabahar and Ahwaz is too large 
with respect to other cities. This high initial investment probably 
has negative impact on the users’ decision to install CCHP systems. 
By considering two techniques this problem may be solved. Firstly, 
thermal storage system beside the solar heating systems or elec¬ 
trical storage system beside the photovoltaic cells can be integrated 
with the proposed CCHP system. This integration can result in 



0 50 100 150 200 250 

Enom(kW) 

- Tehran - Ahwaz 

- Chabahar — — Kamyaran 

Fig. 7. Impact of engine size (E nom ) on the annual S^ 3 for 5 climates. 
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^—Tehran Ahwaz ^—Bandar Anzali ^—Chabahar ^—Kamyaran 
Fig. 10. Impact of engine size (£ n0 m) on the annual S^ 3 for 5 climates. 
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Table 6 

The corresponding Ej{ om for each sub-criterion in the 5 climates 


Conditions 

Elm (kW) 




Tehran 

Ahwaz 

Bandar 

Anzali 

Chabahar Kamyaran 

Maximizing of FESR 

58 

130 

58 

104 

58 

Maximizing exergy 

THTG a 

THTG 

THTG 

THTG 

THTG 

for NPV>0 

63 

165 

66 

148 

59 

IRR > r 

63 

165 

66 

148 

59 

PB < 6 years 

82 

217 

86 

206 

77 

Maximizing NO x reduction 

137 

366 

155 

236 

140 

Maximizing C0 2 reduction 

58 

115 

58 

93 

57 

Maximizing CO reduction 

102 

266 

99 

174 

64 


3 The higher the engine size the greater the parameter, it is calculated according 
to the paragraph after Eq. (33C). 


Table 7 

The optimum engine size, benefits, initial investment and boiler size for 5 climates. 


Parameters 

Tehran 

Ahwaz 

Bandar 

Anzali 

Chabahar 

Kamyaran 

Err m (kw) 

83.00 

211.00 

88.00 

161.00 

78.00 

FESR (%) 

36.23 

35.27 

34.69 

40.08 

37.38 

EXIR (%) 

90.24 

95.57 

90.32 

95.47 

90.05 

NPV (10 7 Rials 3 ) 

14.33 

30.04 

16.12 

7.12 

13.98 

IRR (%) 

16.60 

16.37 

16.67 

15.44 

16.61 

PB (years) 

5.96 

6.04 

5.94 

6.39 

5.96 

NO x Reduction ratio (%) 

99.77 

99.77 

99.77 

99.77 

99.77 

C0 2 Reduction ratio (%) 

45.79 

43.51 

44.71 

46.98 

46.76 

CO Reduction ratio (%) 

80.21 

80.11 

80.19 

80.25 

80.31 

I 0 (10 7 Rials) 

142.09 

346.97 

152.72 

257.90 

137.27 

PVP = 100 x NPV/Io 

10.09 

8.66 

10.56 

2.76 

10.18 

Auxiliary boiler 

73.74 

212.50 

91.62 

102.33 

85.44 

capacity (kW) 







3 10 7 Rials = 1 Million Toman (MT) and officially 1 USD = 12260Rials. 


decreasing the engine and auxiliary boiler sizes while conserving or 
magnifying the other benefits of the CCHP system and as a result 
the initial investment will decrease. These evaluations need more 
calculations and code generation, and are going to be done by the 
authors of this article. Secondly, the government benefits from 
using the CCHP systems, because fuel consumption is lowered, air 
pollution is decreased, these means saving considerable costs for 
the government, therefore government should support this tech¬ 
nology economically by giving subsides and long payback period 
loans to the customers to encourage them for using CCHP systems. 

The cooling system of CCHP system has considerable impact on 
the MCSF. In this research the aim is to propose a sizing technique 
for the prime mover. Therefore, simply the cooling system is 
assumed to be an absorption chiller. However, for comprehensive 
optimization of equipment selection and sizing of cooling system 
there are four main scenarios as follow: ^To use a big absorption 



Tehran Ahwaz Bandar Anzali Chabahar — — Kamyaran 


Fig. 11. Impact of engine size (E nom ) on the auxiliary boiler capacity for 5 climates. 
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Tehran Ahwaz Bandar Anzali Chabahar — — Kamyaran 
Fig. 12. Impact of engine size (£ n0 m) on the initial investment of CCHP for 5 climates. 

chiller with an auxiliary boiler, or direct fire (if the recoverable heat 
of engine is not enough) to provide the total cooling load by a single 
absorption chiller. ^To use two absorption chillers. The first 
operates by the recovered heat from the engine and the second 
operates by an auxiliary boiler or direct fire for the peak time of 
cooling. ►To use a combination of absorption and electrical 
chillers, the absorption chiller uses the recovered heat and the 
electrical chiller uses the electricity of the engine or grid. ►To use 
only electrical chiller to provide the total cooling load. 

The best option depends on the following parameters: Heating, 
cooling and electrical demand of the building, the electricity and 
natural gas tariffs, the economical criteria, thermodynamical eval¬ 
uations and environmental evaluations. This evaluation is going to 
be done by the authors of this paper. 

4.5. Strategy selection results 

As discussed previously in the literature the classical strategies 
of FTL and FEL are usually used as the reference to size the prime 
mover. This method is weak, because only two engine sizes are 
tested to see if which one is better from different points of views. 
This research introduces three more strategies of LFEL, FEPTL, and 
UFTL (Table 3). According to the analyses of this research, strategy 
should not be considered as a constraint. Strategy is a byproduct of 
sizing technique. Fig. 13 shows a comparison between engine size 
proposed by FEL, FTL and MCSF. It shows that MCSF recommends 
using FEPTL for all of the climates. In this strategy, attention is paid 
to maximizing all of the thermodynamical and environmental sub¬ 
criteria. In addition, the project is economically designed to be 
profitable, safe, and have a suitable payback period. Using FTL or 
UFTL increases the I 0 and pollution, and decreases the FESR. In 
addition, using LFEL or FEL is not economical, increases pollution 
and decreases the FESR and EXIR. Therefore, as a result of sizing 
technique (MCSF) the best strategy is FEPTL. It means that the 
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Fig. 13. Strategy proposed by the MCSF versus FTL and FEL for 5 climates. 
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electrical load is fully provided by the CCHP, while the thermal load 
is partially provided and using an auxiliary boiler is necessary at 
this stage. Using auxiliary boiler can be canceled by using solar 
heating or heat storage systems. But using these technologies need 
more attention and calculation which is of interest of authors and is 
going to be done in the future. 

5. Conclusions and further outlooks 

A multi criteria sizing method is introduced to size the prime 
mover of a CCHP system. It uses the FESR and the EXIR as the 
thermodynamical sub-criteria, the NPV, the 1RR and the PB as the 
economical sub-criteria, and the NO x , CO and CO 2 reduction ratio as 
the environmental sub-criteria. Also by using the AHP method, a 
weight is given to every criterion and sub-criterion with respect to 
others. The MCSF is used to size the prime mover of a residential- 
CCHP in five climates. The analyses show that this method is able 
to determine the engine size while satisfying the thermodynamical, 
economical, and environmental constraints and conditions. In 
addition, this technique is able to introduce the best strategy for 
every climate, while the classical methods of sizing do not have this 
ability. The MCSF is easily programmable and user-friendly. The 
CCHP proposed here can be equipped with a thermal storage sys¬ 
tem and an auxiliary solar water heater to decrease the engine and 
auxiliary boiler sizes, increase FESR and possibly decrease of Io. In 
addition using technologies such as photovoltaic cells and electrical 
storage system can magnify the profits of the CCHP system and 
decrease the engine size. According the results, the Io for some 
climates is high, in these cases the government economical sup¬ 
ports is recommended to encourage customers to use CCHP units. 
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